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We consider a bipartite quantum object, composed of a quantum system and a quantum actuator 
which is periodically reset. We show that the reduced dynamics of the system approaches unitarity 
as the reset frequency of the actuator is increased. This phenomenon arises because quantum 
systems interacting for a short time can impact each other faster than they can become significantly 
entangled. In the high reset-frequency limit, the effective Hamiltonian describing the system’s 
unitary evolution depends on the state to which the actuator is reset. This makes it possible 
to indirectly implement a continuous family of effective Hamiltonians on one part of a bipartite 
quantum object, thereby reducing the problem of indirect control (via a quantum actuator) to the 
well-studied one of direct quantum control. 


I. INTRODUCTION 

Coherent control of quantum systems is a central in¬ 
gredient for most implementations of quantum comput¬ 
ing and metrology. In practice, however, interactions be¬ 
tween a controlled quantum system and its environment 
typically lead to a gradual loss of coherence in the system. 
It is, therefore, desirable to perform control operations 
rapidly compared to the system’s coherence times, so as 
to minimize environment-induced errors. This goal poses 
a dilemma, since by coupling a system more strongly to 
a classical controller in order to steer it more quickly, the 
system may also become coupled more strongly to the 
general environment, which can then increase its rate of 
decoherence. Conversely, a system with longer coher¬ 
ence times typically interacts more weakly with all of its 
surroundings, including a classical controller, making it 
difficult to control rapidly. 

An appealing compromise is to employ indirect con¬ 
trol with a setup involving two quantum systems: the 
to-be-controlled quantum system S which is relatively 
well isolated, coupled to a quantum actuator A which is 
classically controlled and interacts more strongly with the 
outside world. This strategy has attracted considerable 
attention since it was found that control of A can yield 
universal control of S [1, 2]. Indirect control schemes 
have shown promise in a variety of settings, including 
spin chains [3-5], superconducting qubits [6], nanome¬ 
chanical resonators [7, 8], and perhaps most notably, in 
nuclear/electron spin systems [9-14], where they are com¬ 
monly used. There is, however, no general recipe for 
mapping a desired unitary on S to a series of control op¬ 
erations on A which will produce it [7]. This is largely 
due to the entangling nature of the system-actuator cou¬ 
pling; a generic operation on A leaves it entangled with 
S , and thus has a net non-unitary effect on the latter. 

Physical settings in which a nuclear spin ( S ) couples 
to an electron spin (A), such as nitrogen-vacancy (NV) 
centers and electron spin resonance (ESR) systems, of¬ 
fer an exceptionally convenient way of surmounting this 
difficulty, as their Hamiltonians can often be cast in the 


form 
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where {|j)} are orthogonal actuator states. The spe¬ 
cial structure of (1) allows one to evolve S unitarily by 
any H$\ conditioned on the state of A [11-14]. This 
approach—which sidesteps the issue of system-actuator 
entanglement by working only with states {|J)} of the 
latter—reduces the problem of indirect quantum control 
to the well-studied one of direct control; i.e., of synthe¬ 
sising a desired unitary using a set of available system 
Hamiltonians [15-20]. 

In this paper, we present an explicit scheme for indi¬ 
rect quantum control which employs a resettable actuator 
to produce unitary evolution of the system, conditioned 
on the actuator’s state. Crucially, our scheme is com¬ 
pletely general in that it can be applied for any S and A. 
Furthermore, our scheme does not require any particular 
form of the <S-A Hamiltonian (e.g., we are not restricted 
to Hamiltonians of the form (1)), but rather, relies on 
the ability to reset A rapidly. In this sense, it reduces 
the general problem of indirect control—for any hybrid 
system—to the much simpler one of direct quantum con¬ 
trol. 


II. THE SCHEME 

Consider a quantum system S coupled to a quantum 
actuator A. Suppose the pair is initially in the state 
p(0) = ps (0) (g) pa, and that A is periodically reset to 
PA at time intervals 5t. (Equivalently, one could couple 
a succession of fresh actuators to S, each prepared in the 
state pa-) For convenience, we will assume the resets to 
be instantaneous, although we shall discuss more realistic 
resetting later. Between successive resets, <S-A evolves as 

p{r) = L{t) P (t), ( 2 ) 

where the superoperator C{t) is known as a Liouvillian. 
We decompose C(r) into parts describing the free dy- 
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namics of S and A, and the coupling between them: 

£(t) = £s + C A + g(,T/St)C SA - (3) 

Here £ 5 , £, 4 , and CsAi which act non-trivially on 
S, A, and <S->4 respectively, are assumed to be time- 
independent and of Lindblad type [21], We include 
an arbitrary piecewise-continuous switching function g : 
[0,1] —y K. in (3) to demonstrate how time dependence 
can be incorporated into this scheme. To facilitate com¬ 
parison with existing indirect control techniques, we will 
be particularly interested in the case where <S->4 evolves 
unitarily between resets, i.e., where £(r) • = •], 

and H(t) = H s + H A + g(T/St)H SA . 

We encode the effect on S of each evolve-and-reset cy¬ 
cle of A in the dynamical map $(<!)£), which acts on the 
system’s density operator as 

rSt 

I dr£(r) 

0 

where T denotes the time-ordering operator. We may 
then write the system’s state after the first cycle as 
ps {St) = $>{5t)ps( 0). If A is reset n times in the interval 
[0, t] (so that t = nSt, where n £ N), the reduced state of 
<S at time t is ps(t) = Q(t/n) n ps(0), where 4>(f/n)" rep¬ 
resents n successive applications of the channel <F(£/n). 

Our scheme utilizes cycles which are short as compared 
to the natural dynamics of <S-M. We will therefore seek to 
expand ps{t) in powers of 1 /n, a small number when the 
actuator is reset at a high rate. In terms of the dynamical 
map <&(t/n) n , we wish to find a series of the form 

$(t/n) n = Ho(£) -I— f£(f) H —2 ^ 2 (t) + ■ • •, (5) 

n n z 


( • ® PA) >, (4) 


$(££) • = tr _4 < T exp 


The key observation is this: in the limit, only l£(0) 
contributes while any 0(59 2 ) terms, in contrast, are sup¬ 
pressed as n —> 00 [22, 23]. 


III. MAIN RESULT 

Remarkably, the fi 0 (t) term in Eq. (5) can be evalu¬ 
ated analogously for any system coupled to a quantum 
actuator that is repeatedly reset. Specifically, Chernoff’s 
theorem [24] (p. 241, see also [25, 26]) gives 

H 0 (£) = lim $(£/n) n = e $lt . ( 8 ) 

n—>o o 

This theorem requires that $ be a continuous function 
of linear contractions on a Banach space, with $(0) = I. 
We verify in the Supplemental Material [27] that $, as 
constructed, satisfies these conditions under the induced 
trace norm on the set of self-adjoint trace-class operators 
acting on the system’s Hilbert space. 

From Eq. ( 6 ), we have that 

$1 • = £5 • +9 tm [Csa{ • ® Pa)} , ( 9 ) 

where g = Jq g(() d( gives the average coupling strength 
between S and A. Observe that to leading order in 1/n, 
the system evolves as ps(t) = e <s>lt ps( 0 ); therefore, $1 
represents an effective Liouvillian for S in the limit of 
frequent actuator resets. We now focus on the case where 
S-M is closed (i.e., where the bipartite object evolves uni¬ 
tarily between resets of A). Eq. (9) simplifies to 

$1 . = -i[i7 eff , •], ( 10 ) 


where each f lk(t) is a superoperator that does not de¬ 
pend on n. Noting that C(5t-() = Cs + C A + g{()Cs A is 
independent of St (since the argument of g is scaled by St, 
see Eq. (3)), we proceed by first expanding <j>(5f) asymp¬ 
totically for small St as $>{5t) = I + St <I>i + St 2 $ 2 + ..., 
where 


K 1 


r(k-i 


• = tr _4 / d<i 


dC2 • • 


dCfc 


K 

n - Cj) (• ®Pa) ( 6 ) 


m 


is a superoperator with no St dependence. 

Before presenting our main result, let us consider an 
analogous but simpler situation, namely the case of a 
matrix R{9) £ SO(2) representing a rotation by 9 in 
R 2 . The action of R{9) can be obtained as the out¬ 
come of a series of infinitesimal rotations, each given by 
I + S9 R'(0) + O{50 2 ), for an arbitrary 0(S9 2 ) term. Con¬ 
cretely, setting 59 = 9/n: 


lim 

n—>o o 


I+-R'(0) + O 

n 



= e R'(o)0 


R{9). (7) 


where 


H eS = H s + g tr _4 (H sa Pa) ■ (H) 

We show explicitly the details of this calculation in the 
Supplemental Material [27]. Clearly H' ef{ = H e g; it fol¬ 
lows that, to leading order, the system’s reduced dynam¬ 
ics is unitary, and described by 

= [H e s, ps{t )]■ ( 12 ) 

If, for example, the interaction Hamiltonian has the 
form Hsa = B (g) C, the system’s reduced dynamics 
will be well-described by H e g = Hs + g(C)B, where 
(C) = tr{CpA), when the reset rate is high. More gen¬ 
erally, different types of coupling between S and A will 
lead to different effective Hamiltonians, generating uni¬ 
tary system dynamics conditioned upon p A - Therefore, 
by gradually varying the state in which A is prepared, 
one can use this scheme to implement entire families of 
Hamiltonians on S , reducing the problem of indirect con¬ 
trol to one of direct quantum control. Moreover, for 
finite-dimensional systems, our scheme gives universal 
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unitary control of S provided H$ and H$a are not re¬ 
lated by some symmetry [15]. 

One can view S as receiving a discrete kick from A 
with each cycle. Notice that &(5t) —> I as St —>• 0; in 
other words, as the reset rate becomes large, the kicks to 
S become weaker but more frequent. These competing 
trends underpin our scheme: When the cycles are short, 
<3>(<$t) ss I + St $i, i.e., the impact of each kick is mostly 
encoded in <f>i. As in (7), the aggregate effect of many 
such kicks depends only on d>i to leading order. Cru¬ 
cially though, d>! generically describes a non-trivial, but 
non-entangling operation on <S-A. This is not unique 
to our scheme: the non-entangling nature of the first- 
order Dyson series term has been used in very different 
contexts, see, e.g., [28]. The present scheme, therefore, 
exploits the general phenomenon that two quantum ob¬ 
jects interacting for a short time St can impact each other 
faster than they can become significantly entangled. By 
mimicking a series of short interactions between S and A 
through frequent resets of the latter, we can effectively 
modify the system’s Hamiltonian without significantly 
increasing the system’s entropy. Along the same lines, 
one can also think of the frequent actuator resets as serv¬ 
ing to keep A in its initial state. The resulting effect 
is for S to evolve according to the portion of the total 
Hamiltonian acting on the system’s Hilbert space, up to 
corrections of order 0 (l/n). 

IV. EXAMPLE 

To illustrate our scheme, we consider an archetypal 
problem in quantum control: that of steering a harmonic 
oscillator S indirectly via a d-level actuator A. Exist¬ 
ing techniques depend crucially on d and/or on the na¬ 
ture of the <S-A coupling. For instance, [29] requires a 
Jaynes-Cummings (JC) interaction with d = 2, [30] re¬ 
quires a JC-like coupling with d = 3, and [7] requires 
switching between two distinct interaction Hamiltonians. 
Our scheme, in contrast, is model-independent. For illus¬ 
tration, however, we pick d = 2 and an 5-A Hamiltonian 

H(t) = Hva)a + oz + hg{r/ 8 t)X <g> n ■ er, (13) 

where X = [a + a^)/ 2 is a quadrature operator on S and 
cr = ( a x ,ay,cr z ). For different choices of n, Eq. (13) de¬ 
scribes several hybrid objects amenable to indirect quan¬ 
tum control; for example, a nanomechanical resonator 
coupled to a superconducting qubit [7, 31, 32] or to an 
electron spin [ 8 ], as well as circuit (cavity) quantum elec¬ 
trodynamics setups, where an electromagnetic field mode 
couples to superconducting qubits (atoms). 

A direct application of Eq. (11) yields an effective 
Hamiltonian 

H c g = hva)a + hg{\\ • er)X, (14) 

on S , which is manifestly dependent on pa- It follows 
from Lloyd’s well-known argument in Ref. [15] (see also 


[33]) that any Hamiltonian of the form 

H = cia 1 a + C 2 X + C 3 P, (15) 

can be implemented on S by suitably varying the state 
to which A is reset, where P = (a — a^)/2 i. 

In implementations, the accuracy with which one can 
enact evolution by Lf e ff through frequent resets of A will 
be important. To quantify this accuracy, we consider an 
initial system state /?s(0) and evolve it numerically for a 
time t according to the full Hamiltonian (13), where A 
is reset at a rate f = n/t. We then compute the fidelity 
F (see caption of Fig. 1 for definition) as a function of 
t, between the resulting reduced system state and the 
time-evolved state e~ ltHBi! ^ h ps( 0 )e ltHB!! ^ h that we wish 
to obtain. The resulting fidelity is shown in Fig. 1 for 
three different reset rates. 



FIG. 1. (Color online) Deviation between 77 e g and full S-A 
evolution. Shown is the fidelity F(t) = yj | ps(t) \ ip(t)) 

between \ip{t)) = exp(—itH e s/h)\a) and ps(t), the reduced 
state of S under evolution by (13), for three different resetting 
rates /. (From top to bottom: / = 10u, / = 5n, f = 2v.) 
The initial state of the system ps( 0) = |n)(a| is coherent 
with a = (1 + i)/\/2. We set v = u, n = (1,0,0) and pA = 
(1 + ct t )/2 for illustration, and chose the switching function 
as g(r/St) = 2u sin 2 (nr/St) [34]. 

We observe that two qualitatively distinct types of er¬ 
ror arise in implementing H e g on S with our scheme, 
see Fig. 1: (i) slowly growing deviation from F = 1 as t 
increases (apparent in the main panel), and (ii) fast “wig¬ 
gles” on these otherwise slowly varying curves (shown in 
the zoomed inset). Type (i) error, which we call dissipa¬ 
tive error , arises because S and A get slightly entangled 
during cycles of finite duration. Over time, this entan¬ 
glement causes information about S to become lost due 
to the actuator resets, contributing a non-unitary compo¬ 
nent to system’s evolution, which accumulates with every 
cycle. Thus, dissipative error becomes important on long 
timescales. 

Type (ii) error, which we name stroboscopic error , also 
arises from a finite reset rate /; however, it is important 
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only on comparatively short timescales. When / < oo, 
our scheme approximates smooth evolution by H e g using 
discrete non-unitary kicks from A, which can have com¬ 
plex effects on S. Our asymptotic analysis has been fo¬ 
cused on the total effect of each kick, described by Eq. (5) 
for integer n. However, the system can display compli¬ 
cated non-unitary dynamics during a cycle of finite dura¬ 
tion, the details of which are not described by H e g (nor, 
more generally, by Qi(t), 0 2 (f), ...). Stroboscopic error, 
then, is the temporary deviation during each cycle be¬ 
tween the full open dynamics of S and the smooth path 
described in (5). Because this type of error vanishes at 
the end of every cycle, it does not accumulate with t, 
and hence, is primarily important over short timescales. 
We will now establish how both the dissipative and the 
stroboscopic error decay as the cycle length goes to zero. 

V. ERROR ANALYSIS 

Dissipative error in implementing H c g arises from the 
fii(f), fl 2 (f), ... terms in (5), which generically introduce 
dissipation into the system’s evolution. These superop¬ 
erators can be expressed as functions of $i, <h 2 , ... by 
adapting the method developed in [35] (Section 4). We 
will use this method to find the form of fii (t), the leading- 
order source of dissipative error. 

Consider the function 

u(t) = $(rt/n) n n 0 [(l - r)t], (16) 

chosen so that v(l) = §(t/n) n and u(0) = O 0 (t). Observe 
that with this definition of v we have 

<E>(f/n)™ — B 0 (f) = u(l) — i>(0) = f w , (r)dr. (17) 

Jo 

One rapidly arrives at expressions for Oi(f), 0 2 (t), ... 
by expanding v'{t) asymptotically in 1 /n and matching 
powers with Eq. (5). In particular: 

-Sh(t) = - [ e* lrt ($ 2 -s $ i) e * l(1 ~ T)t dr. (18) 

n n J o V 2 / 

In terms of the reset rate / = n/t, Eq. (18) scales as 
0(t/f). (This statement is readily formalized by noting 
that the induced trace norm of the integral in (18) is inde¬ 
pendent of t.) The t/f scaling is apparent in Fig. 1, where 
the reset rates considered are sufficiently high that Qi(t) 
is the dominant source of dissipative error. In particular, 
observe that for each / plotted, the deviation between 
the reduced dynamics and the Ff e ff-generated trajectory 
(corresponding to Dg(t), the leading-order term in 1 /n) 
is nearly linear in t. Furthermore, the slope of the lines 
scale inversely with /. 

We now turn our attention to the stroboscopic error. 
The dynamical map in Eq. (5) describes the system’s 
evolution for a time t = nSt, corresponding to an integer 
number of cycles. In other words, it gives ps at the end 


of each cycle. However, between successive resets (i.e., 
mid-kick), S will temporarily stray from the smooth tra¬ 
jectory given in (5). 

One can quantify the stroboscopic error by comparing 
the desired evolution by H e g with the system’s reduced 
dynamics between successive actuator resets. We use a 
truncated Dyson series to perform this comparison: If 
ps(t n ) is the system’s state after cycle n, then during 
cycle n + 1, evolution by H e s would give 

Ps S \tn + t) °= Ps(t n ) - Ps(t n )], (19) 

for 0 < t < 6 t. When the cycles are short, i.e., when 
t < St <C 1 in units set by the largest characteristic 
frequency of H(t ), the 0(t 2 ) terms in the Dyson series 
will be subdominant, and so we work only to first order 
in t. 

Eq. (19) is the evolution our scheme seeks to imple¬ 
ment. However, the actual full evolution of S between 
successive resets is computed by evolving <S-_4 according 
to H(t ) and then tracing out the actuator: 

Ps^ttn+T) °= ) ps(tn)~^tl A [ J H(t')At' , ps (tn) Pa\ ■ 

m 

Concretely, stroboscopic error is described by the differ¬ 
ence between Eqs. (19) and (20): 

P S ( tn+T)-p y s '{t n + T ) = (21) 

St f T 1 

9~ — J Q 5(0 dC \ [tr A (H SA p A ), ps{t n )\. 

The braced term in (21) is bounded above in magnitude 
by 2<7 max (<5t — t)/t, where (? m ax is the largest coupling 
strength attained in each cycle. (Thus, stroboscopic er¬ 
ror is reduced when the <S->4 coupling remains weak.) 
The commutator, in contrast, can vary arbitrarily with 
t n , depending on the nature of H(t). However, it is al¬ 
ways suppressed by a prefactor of St — r, and so Eq. (21) 
generically scales as 0(St — r), which is upper-bounded 
by 0(dt). In terms of reset frequency then, stroboscopic 
error scales as 0 (1/f). 

Unlike dissipative error, which accumulates with t, 
stroboscopic error vanishes with each reset. Thus, while 
the former type is reduced by implementing H e Q for short 
durations t, the latter can be sidestepped entirely by 
choosing / and t to give an integer number of cycles. In 
general, both types of error can be suppressed arbitrarily 
by choosing an / that is large compared to characteristic 
frequencies of S and A. 

VI. DISCUSSION & OUTLOOK 

Our scheme is somewhat reminiscent of the quantum 
Zeno effect (QZE) [36-38]; both our scheme and the QZE 
feature emergent unitary dynamics generated by an ef¬ 
fective Hamiltonian, which results from rapidly repeated 
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operations. However, our scheme does not invoke any 
type of measurement, nor any notion of “leakage” be¬ 
tween measurement eigenstates. Moreover, the QZE typ¬ 
ically involves repeated kicks which remain strong (or 
at least non-vanishing) in the high-frequency limit [39]. 
In contrast, the effect on S of each individual cycle in 
our scheme vanishes as / —> oo. Zanardi and Cam¬ 
pos Venuti recently discovered a phenomenon closely re¬ 
lated to the QZE, wherein unitary evolution, generated 
by a “dissipation-projected Hamiltonian”, can arise in an 
open system. While their results and ours may both stem 
from a common fundamental principle, the scaling that 
underpins Refs. [40, 41] is entirely absent in our scheme. 

While we have assumed for illustration that A is reset 
instantaneously and at regular intervals, neither of these 
idealizations are essential to our scheme. Notice from 
Eq. (9) that does not contribute to $i; therefore, 
S would still evolve unitarily by H c g to leading order 
if A were reset through a gradual non-unitary process. 
We note also that Chernoff’s theorem —which gives the 
evolution of S in the high / regime—can be generalized 
to describe cycles of non-uniform duration. Specifically, 
when n actuator resets are performed in a time t, the 


system’s evolution will be well-described by H e s (or more 
generally, by $i) when the longest time between actuator 
resets is sufficiently short [42, 43]. 

Finally, we wish to compare our control scheme with 
the one proposed in [1]: Whereas we employ a resettable 
quantum actuator to achieve indirect unitary control of a 
system, Lloyd and Viola used a resettable and measurable 
ancilla to achieve arbitrary open-system dynamics. The 
similar requirements of both schemes suggest the possi¬ 
bility of implementing arbitrary open-system dynamics 
on S through indirect control. 
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